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A  Diffusion  Defined  on  A  Fractal  State  Space 


William  Bernard  Krebs 

Abstract:  We  define  a  fractal  in  the  plane  known  as  the  Vicsek  Snowflake  by  constructing  a 
skeletal  lattice  graph  and  then  rescaling  spatial  dimensions  to  give  a  sequence  of  lattices  that 
converges  to  a  fractal.  By  defining  a  simple  random  walk  on  the  skeletal  lattice  and  then 
rescaling  both  time  and  space,  we  define  a  sequence  of  random  walks  on  the  approximating 
lattices  that  converge  weakly  to  a  limiting  process  on  the  snowflake.  We-show-th&t^this 
limit  has  continuous  sample  paths  and  the  strong  Markov  property,  and  that  it  is  the 
unique  diffusion  limit  of  random  walk  on  the  snowflake  in  a  natural  sense.  We-show-thatf — 
this  diffusion  has  a  scaling  property  reminiscent  of  Brownian  motion,  and  we  introduce 
a  coupling  argument  to  show  that  the  diffusion  has  transition  densities  with  respect  to 
Hausdorff  measure  on  the  snowflake.  , 


Keywords:  Diffusions,  fractals,  <  ■.  (•  ts  ' — 
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1.  Introduction: 

Construct  a  figure  in  the  unit  square  bv  the  following  recursive  procedure.  Let  t/o  denote  the  unit  square. 
Construct  Gi  by  deleting  from  Go  four  squares,  each  with  edge  length  |,  centered  along  the  four  edges  of  Go- 
G 1  will  consist  of  five  squares  with  edge  |  whose  corners  overlap.  At  stage  n,  Gn  will  consist  of  5n  squares, 
each  with  edge  length  3“".  To  construct  Gn+i  from  Gn,  take  each  square  S  composing  G„ ,  and  delete  the 
four  square  centered  along  the  edges  of  with  edges  of  length  3-n_I.  Gn+ 1  then  consists  of  5n+1  squares 
with  edges  of  length  3_n_1. 
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Take  G  =  Some  easy  topology  shows  that  G  is  a  closed  connected  set.,  with  Lebesgue  measure 

0.  In  fact,  if  is  not  hard  to  show  that  G  has  finite  Hausdorff  log3  5-dimensional  measure.  In  the  spirit  of 
Mandelbrot,  G  is  a  fractal  with  starter  polygon  Gi-  Extensive  treatments  of  such  fractal  sets  have  been  gi\en 
by  various  authors.  (See,  for  example  Hutchinson[12]  or  Barnsley  and  Demko[4]). 

A  number  of  authors  have  treated  the  problem  of  constructing  diffusions  on  nested  fractals.  Particular 
attention  has  been  paid  to  diffusions  on  the  Sierpinski  gasket,  a  fractal  constructed  from  a  unit  equilateral 
triangle  by  successsively  deleting  “middle”  triangles.  Goldstein[10]  and  Kusuoka[13]  constructed  a  Brownian 
motion  on  the  Sierpinski  gasket,  using  a  decimation-invariance  property.  Barlow  and  Perkins[3]  have  studied 
this  Brownian  motion  comprehensively.  Browniand  motion  on  the  Sierpinski  gasket  is  broadly  similar  to  the 
natural  diffusion  on  the  Vicsek  snowflake,  and  the  results  of  these  authors  are  generally  similar  to  those  in  the 
present  work.  I  fully  acknowledge  the  priority  of  their  results.  More  recently,  Lindstrom[14]  has  constructed 
a  Brownian  motion  on  any  fractal  set  satisfying  a  general  set  of  nesting  axioms  from  a  sequence  of  random 
walks,  provided  that  the  distribution  of  the  random  walk  satisfies  a  non-degeneracy  condition. 


The  first  objective  of  this  paper  is  to  construct  a  diffuion  on  the  Vicsek  snowflake,  starting  from  an  non¬ 
degenerate  random  walk  model.  In  an  important  respect,  the  problem  of  defining  diffusions  on  the  snow  flake 
is  more  complicated  than  defining  diffusions  on  the  Sierpinski  gasket.  On  the  snowflake,  one  can  define  a 
variety  of  random  walk  models  that  are  symmetric  under  the  natural  symmetries  of  the  square.  A  natural 
question  is  whether  one  can  constiuct  a  diffusion  for  any  such  model.  Another  is  whether  the  diffusion  on 
the  fractal  is  unique  such  diffusion  are  unique.  The  snowflake  seems  to  be  the  simplest  nested  fractal  where 
such  questions  arise.  For  the  snowflake,  the  answer  is  that  if  the  random  walk  is  not  degenerate  then  the  a 


unique  diffusion  limit  exists  independent  of  the  underlying  random  walk  model.  The 
for  general  nested  fractals  remains  unsolved  at  the  time  of  this  writing. 
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2.  Constructing  the  Diffusion: 

Consider  the  following  system  of  transformations: 

Mi  :  x  — f  3-1z  Mu  :  x  — i  (3-1  ■  x)  +  (2,0) 

Mi  -.x-tZ-1  .*  +  (1,1)  — (3_1-*)  +  (0>2)  (2.1) 

M3  :  x  -f  3"1  •  x  +  (2, 2) 

By  inspection,  Mi}  ..,M5  are  strict  contractions,  with  fixed  points  (0,0),  (l^,l|),  (3,3),  (3,0)  and  (0,3), 
respectively.  For  bounded  subsets  A  of  R2,  define  M(A)  =  u\M,(A).  It  is  well-known  to  geometers  that  the 
transformation  M  has  a  unique  compact  invariant  set,  whose  Hausdorff  dimension  may  easilj  be  computed 
as  log35.(  See,  for  example,  Barnsley  and  Demko  [4],  Dubins  and  Freedman  [6]  or  Hutchinson  [12]  ).  We 
will  call  this  set  the  bounded  Vicsek  snowflake ,  and  denote  it  by  Tj. 

For  future  reference,  we  establish  the  following  definition.  Let  M,, , ...,  M,n  be  some  sequence  of  the  trans¬ 
formations  M\  ...,M5.  Let  S  =  M,-,  o  ...  o  M,n(r»).  We  will  call  S  a  square  of  r^.  We  will  also  need  an 
unbounded  version  of  our  state  space.  Let  T  =  U£L03nF&.  T  also  has  Hausdorff  dimension  log3o,  and  has 
the  property  that  =  T.  We  will  call  F  the  unbounded  Vicsek  snowflake. 

We  wish  to  construct  a  diffusion  process  on  T  and  study  its  basic  properties.  We  shall  do  this  b\  defining 
random  walks  on  a  sequence  of  lattices  that  approximate  T,  which  we  new  construct. 

In  the  unit  square,  let  U  denote  the  complete  graph  on  the  corners  of  [0,  l]2.  Let  V(LT)  and  £(U)  denote  the 
vertices  and  edges  of  U,  respectively. 

Define  a  new  graph  U0  with  vertex  set  V(U0)  and  edge  set  £{U0)  by  taking 

V(f'o)  =  V(U)  U  (V(L7)  +  (1, 1))  U  (V(t/)  +  (0, 2)]  u  [V(U)  +  (2,0)]  U  (V(t0  +  (2, 2)] 

(2.2) 

S(U0)  =  £{U)  U  [£(U)  +  (1, 1)]  U  [£((/)  +  (0,2))  U  [f(l/)  +  (2,0)]  U  (f(l/)  +  (2,2)] 
where  the  arithmetic  is  done  componentwise.  We  cal!  Uo  the  unit  snouflake  lattice.  Note  that  U0  lies  in  the 

square  [0,3]2;  we  will  call  the  points  (0,0),  (0,3), (3, 3),  and  (0,3)  the  outer  corners  of  Uo. 
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P[A'„+i  =  !/]  =  p  if  x  and  y  are  diagonally  adjacent 

P[Ar„+i  =;/]  =  (  1  —  p)/ 2  if  x  and  y  are  vertically  or  horizontally  adjacent 

If  Nr  =  6: 
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P[ A',,+i  =  y]  =  p/2  if  x  and  y  are  diagonally  adjacent 

P[A'„+i  =  t/]  =  (1  -  p)/4  if  x  and  y  are  vertically  or  horizontally  adjacent 

Here,  0  <  p  <  1  is  a  fixed  but  arbitrary  parameter.  This  defines  a  random  walk  on  G,  which  we  denote  by 
X'P.  Thus,  if  Xg  =  (0,0)  then  P[ X\  =  (1, 1)]  =  p  while  P[X f  =  (1,0)]  =  (1  -  p)/ 2.  If  X[  =  (1, 1),  then 
P[Xl  =  (0, 0)]  =  p/2  and  P\Xg  =  (1, 0)]  =  (1  -  p)/4  =  P[Xg  =  (2, 1)] 

We  begin  by  studying  the  special  case  where  p  =  (1  —p)/2  =  |  ,  and  Xg  =  0.  This  will  define  a  simple 
random  walk  on  the  graph  G,  starting  from  0.  We  call  this  discrete-time  Markov  chain,  which  we  will  call 
random  walk  on  the  snowflake  lattice,.  For  convenience,  we  will  write  X1/3  =  A”,  o 

Let  and  T2  be  the  sequences  of  times  between  visits  by  A'n  to  distinct  points  of  3G  and  3 -  G  respectively. 
Since  32  G  C  3  G,  T3  =  •  %  the  nesting  property  of  the  lattice,  the  distribution  of  b(j)-a(j)  is  the 

same  as  the  distribution  of  T/ ,  and  the  Markov  property  of  A'n  shows  that  b(j)-a(j)  and  T/,k  =  a(j), ... ,b(j ) 
are  independent,  are  independent  random  variables,  equal  in  distribution  to  T.  Thus,  for  each  j  Tf  has  the 
distribution  of  the  second  generation  of  a  branching  process  with  offspring  distribution  equal  to  that,  of  T. 

Similarly,  let  T"  be  the  times  between  visits  to  distinct  points  of  3"  •  G.  A  similar  argument  shows  that,  for 
each  I,  Tj/  has  t’  '  distribution  of  the  nth  generation  of  a  branching  process,  again  with  offspring  distribution 
equal  to  that  of  T. 


Let  /(«)  be  the  generating  function  of  T.  By  direct  calculation,  we  can  show  that 


/(»)  = 


(3  -  2u)(12  -  12u  -f  ji2)  ’ 


(See  Section  4  for  the  details  of  the  computation.)  Differentiating  /  shows  that  ET  =  lo  and  Var(r)  =  114, 
so  the  branching  process  {7’n}  is  obviously  supercritical.  Since  T  is  always  strictly  greater  than  0,  {T")  has 
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extinction  probability  0.  Since  this  branching  property  of  the  hitting  times  of  Xn  plays  a  key  role  in  the 
remainder  of  this  section,  we  will  review  some  standard  theory  of  branching  processes. 

Theorem  2.1.  Let  Z„  be  a  branching  process,  with  Zo  =  1  and  let  f  be  the  generating  function  of  the 
offspring  distribution.  Suppose  that  1  <  f'(l-)  =  m  <  oc  and  /"(!—)  <  oo.  Let  W„  =  Znm~n.  Then 
there  exists  a  random  variable  IV  with  EW  =  1  such  that  lVn  —  IV  a.  s.  and  in  L~  P[ IV  =  0]  = 
P[Z„  =  0  for  some  »].  If  <f>(u)  =  e~uW ,  then  satisfies  Abel’s  functional  equation  f>(u)  =  f 

Proof:  This  is  Theorems  1  and  2,  and  equation  (5)  in  Athreya  and  Ney  [2], 

In  particular,  IV„  — *  IV  in  distribution. 

The  theorem  implies  that  lo~nTn  converges  in  distribution  to  some  random  variable  IV,  with  £IV  =  1. 
As  T"  can  never  be  0,  IV  is  strictly  positive  almost  surely,  and  ^(A)  =  Ee~xw  satisfies  p(A)  =  /(9(15_1A). 

For  m  s=  0.  1,...,  define  a  stochastic  process  on  T  by  Ym(t)  =  3“mA’([15mtl).  (  [x]  denotes  the  greatest 
integer  less  than  or  equal  to  z).  Observe  that  for  each  n,  Y„(t)  is  a  random  walk  on  3-ri  •  G.  Let  D[0,oo]  be 
the  set  of  functions  w  :  R+  — ♦  T  that  are  right  continuous  and  have  left  hand  limits  for  all  i. 

Theorem  2.2.  The  sequence  of  processes  On(0)«>o  ‘s  tight  in  D[0,  cc-].  IfY„i(t)  is  a  subsequence  of  {}'„(/)} 
converging  weakly  to  a  process  Yt,  then  Y,  has  continuous  sample  paths. 

The  theorem  follows  from  an  estimate  of  the  moments  of  the  displacements  ]|V„(t)  -  V„(s)||,  which  we  state 
as  a  lemma. 

Lemma  2.3.  For  n  =  1,2,  ...,0  <  s  <1  <  00,7  >  0, 

F||y„(0  -  y„(«)|P  <  2%/?  •  [3~’n  +  C  •  (1  -  s)M] 
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(2.5) 


where  p  =  log15  3  and  C  is  a  constant  independent  of  n. 


Proof:  Let  O  =  {q  E  Q  ■  q  =  p  ■  lb~n,p,q  €  Z}.  Since  Y„(i )  jumps  only  at.  points  in  Q,  it  will  suffice 
to  estimate  2?|pn(?)  —  ^n(r)IP>  f°r  arbitrary  §,r  6  O.  We  make  the  following  displacement  estimate  for 
random  walk  {Xn}  on  the  snowflake  lattice: 

P[||A'„  -  Xm\\  >  2V2  ■  3*]  <  P[Tk  <n-m]  (2.6) 

To  establish  this  estimate,  we  observe  that  2>/2  •  3*  is  the  diameter  of  two  diagonally  adjoining  squares  in 
the  3*  •  G.  If  ||AT„  —  Xm\\  >  2\/2  ■  3*,  then  between  epochs  m  and  n,  A'*  must  visit  two  distinct  points  of 
3*  •  G. 


To  use  estimate  (6),  we  write 


£||yn(?H'„(r)|P 

<  (2v/2 y  •  3 -v>P  [||Tn(?)  -  Vn(r)||  <  2V2 . 3"” 

OO 

+  (2V2y  ■  £  3,-7P  [2V2 . 3f  <  ||y„(j)  -  y„(r)||  <  2V2  ■  3f+1] 

*=-n 

<  {2Viy  ■  3 -^p  [||y„(?)  -  y„(r)||  <  2V2  ■  3-"] 

+  (2 V5y  •  (1  -  3-*)  •  f;  3 hp  [||y„(9)  -  y„(r)||  >  2V2  •  s'] 

i=-n 

<  (2V2y  •  3 -*nP  [||A'(15n  ■  q)  -  A'(15n  ■  r)||  <  2V2} 


OO 

+  (2\/2)7  ■  (1  —  3-7)  •  ^3(*“")‘>p  |||A’(15"  ■  q)  -  A’(15"  ■  r)||  >  2%/2  -  3*1 
Jt=o 

<(2  V2.y  ■  3->n +  (l-3-f).^{k~nhP[lok -n  <  1 5- •(</-»•)] 


k=0 


Inequality  (8)  comes  from  substituting 

P[ 2v/2  •  3'  <||y„(9)  -  y„(r)||  <  2V2. 3*1]  = 

P[2V2  ■  3'  <  ||yn(9)  -  y„(r)||l  -  p[2\/2  •  3f+1  <  ||yn(9)  -  yn(r)||l 


(") 

(8) 


(9) 

(10) 


into  the  estimate  and  rearranging  the  terms.  (9)  follows  from  the  definition  of  Yn(q),  and  (10)  is  from 
estimate  (6). 


To  estimate  P  [l5fc  •  7*  <  15~*  ■  15"  •  ( q  -  r)]  we  first  estimate  the  Laplace  transforms  of  {7*}.  Let  fa  and 
<$>  be  the  Laplace  transforms  of  T-K  and  IT,  respectively.  As  7*  has  the  distributions  of  a  branching  process 
and  15“  *7*  — *  IT,  it  is  not  hard  to  show  that  <f>k  T  <f> ■  So,  it  will  suffice  to  estimate  <j>. 

To  estimate  ©(15*  •  ( q  -  r)-1),  let  h{u)  —  -  log(<i(t/))  be  the  cumulant.  generating  function  of  IT.  h  satisfies 
the  functional  equation  h(u)  =  log  (/  (exp  (A(15-1  •  u)))),  where  /  is  the  generating  function  of  T  =  T1. 
Let  1  <  s  <  15.  Since  h  is  non-decreasing,  we  have, 

5  W  =  108  (/  (esp('l(15'1  ■ "»> 

=  (logfofnji)  ■|0g  ((3 -2«»))(]25-12dW +  «»)-))  (2'U) 

_  ?  l°g((3  ~  2fau))(  12  -  12^(u)  +  d(u)2)) 

log(^(u)) 

The  second  term  goes  to  0  as  u  — oo.  Since  h  is  monotone  increasing, 

limsup/i($u)//i(u)  <  3,  1  <  s  <  15.  (2.12) 

11  — CO 

This  shows  that  h  is  a  function  of  dominated  variation,  which  implies  the  existence  of  constants  C\  and  C2, 
such  that  C\up  <  h(v)  <  Ciup,  where  p  =  log153.  (See  Feller[S],  de  Haan  and  Stadtmuller[ll]).  Restating 
this  in  terms  of  <?  shows  that  fat)  <  exp(-Cif/’). 

We  now  apply  a  standard  technical  result  on  Laplace  transforms. 

Lemma  2.4.  Let  U  be  the  distribution  function  of  a  random  variable,  let  C(u)  be  its  Laplace  transform, 
and  let  a  >  0.  Then,  for  any  t  >  0,  U(t)  <  e~a  •  fat-1). 


Proof:  This  is  proved  in  the  Corollary  to  Theorem  XIII.5.1  in  Feller  [91. 


Applying  these  two  lemmas  to  our  series  gives 


]r3(*-»hJp[l5*  ,Tk  <  15-*  .  15n  .  (9_  r)  1 

kzz  0  J 

<  c ,'15  ■  E  ■  ( q  -  r)"1)  (2.13) 

k- 0 
oo 

<e 1/15  •  E  3(i-n)^(15t-".(g-r)-1)  (2.14) 

*  =  -CO 

00  00 

<  e1'15  •  [E  3-m7  -f  E  3m7  exp  (-15m'  •<?,(?-  r)"')l  (2.15) 

m= ’  m=0 

<  OO 


To  complete  the  proof  of  the  estimate,  let  j(/j  =  [log15il  +  1  for  i  >  0  and  let  F(t)  =  e1/15  •  3^(0  . 
£~-oo  3‘”r^(15l).  Substitute  inequality  (2.14)  into  the  left-hand  term  in  (2.10)  and  let  m  =  k  -  n,  to  get 


£||*n(g)  -  Tn(r)|P  <  (2\/2)'y  • 

<  (2  V2y  ■ 


3-^4(i_3"^).eJ/13.  E  3m^(15m  ■  (q  -  r)"1) 

mr: -oo 

3-1"  +3V(?-r)(1  _3-7).el/l5  .  ^2  3m^(l5m) 


(2.16) 

(2.17) 


From  (15),  3mT'^(15m)  =  C  <  oo.  Since  j(q  -  r)  <  logJ5(g  -  r)  +  1,  3^(«-r)  •  e1/15  <  (?  -  r)^. 

Subsitituting  these  expressions  into  (2.17)  gives 


E\\Yn(g)-Yn(rW< 


3-in  +  (9  —  r)n  •  (1  —  3”'r)  •  E  3mV(15m) 

m=  — 00 


which  proves  the  lemma. 


(2.1S) 


Using  our  estimate,  we  can  deduce  both  weak  convergence  of  }'„.(/)  and  sample  continuitv  of  1).  Weak 
convergence  of  (ln(t))  follows  from  the  standard  result  on  convergence  of  stochastic  processes. 


Theorem  2.5.  Suppose  that 


(-Y„(/i) . A '„(/*)) (A’(/,) . X(tk)) 
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holds  whenever  t\  <  ...  <  are  points  where  X(t)  is  almost  surely  continuous,  and  that 


£[||A'„(<)  -  A'n(i,)|n|A'„(/2)  -  A«(Oir]  <  (F(h)  -  F(U))2° 

for  t\  <  t  <  <2,  n  >  1,  7  >  0,  a  >  and  F  a  nondecreasing  continuous  function.  Then  Xn  —*o  X. 

Proof:  This  is  Theorem  15.6  in  Billingsley  [5]. 

Proof  of  the  main  theorem:  To  show  that  there  exists  a  weakly  convergent  sequence  {}'„/(<)},  it  suffices 
to  show  that  {!„}  is  tight.  Let  s  <  u  <t  and  consider 

E[\\yn(u)  -  v«(*)ipr»(<)  -  y„(«)ip  3  (2.19) 

If  |s  -  f|  <  15-n  then 

^(iii»(«)  -  y»winii'n(o  -  y„(u)in  =  0  (2.20) 

for  any  s  <u<t,  because  Yn{t)  jumps  only  at  integral  multiples  of  15“n,  and  if  |s  -  i\  <  15~n  then  there 
can  be  at  most  one  such  multiple  in  [s,<].  So,  suppose  that.  |s  -  t|  >  15~n.  Apply  Holder’s  inequality  and 
the  monotonicity  of  £||yn(u)  -  yn(s)||T  to  see  that 

^[i|y»(«)-y»wrn(o-y»(«)ip] 

<E\\Yn(t)-Yn(s)\\^ 

<  (2\/2)3^  [3-n2’  +  (1  -  ■  (1  -  s)-^)  (2.21) 

But,  since  |s  -  f|  >  lo"",  3-"2^  <  (s  - 1)-^.  Substitute  this  into  (2.21)  to  get 

£(llyn(«)  -  yn(OIPI|y»(0  -  Vp(«)|P  3  <  2  •  (2%/2)2-»  ■  (/  -  sf”  (2.22) 

This  suffices  to  prove  tightness. 
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Choose  a  subsequence  n'  so  that  j  Yn>(t)  —*d  Yt  ■  To  show  that  the  sample  paths  of  Yt  are  continuous  with 
probability  1,  let  n  — ►  oo  in  estimate  (2.6),  giving 


B\\Yt  -  mP  <<  C  ■  (t  -  s)” 


(2.23) 


Kolmogorov’s  criterion  for  sample  path  continuity  applies,  which  completes  the  proof  of  our  main  theorem. 


Although  we  have  written  our  proofs  for  the  special  case  where  Vo  =  0,  the  proof  will  also  work  for  Yq  =  x 
for  an  arbitrary  x  €  Too,  with. a  minor  modification  of  the  sequences  Y„(t). 


For  n  =  1,2,...  and  m  <  n  let  be  the  time  of  the  kth  visit  of  the  random  walk  Yn(t)  to  Gm  and  let 
T™k  =  S™rnk  —  S™k_1  be  the  kth  interarrival  time  for  S™|Jt.  Previously,  we  observed  that  Yn(S™k)k  =  1,2, ... 
is  a  random  walk  on  Gm,  for  all  m  <  n.  We  can  extend  this  property  to  Yt. 


Proposition  2.6.  Let  Yn>  — ►  Y.  Tn><m< *  converges  in  distribution  to  a  random  variable  T„hk  for  every 
k  >  0;  furthermore,  the  sequence  is  independent  and  identically  distributed  for  all  m.  Yn'{Sn\m,k) 

converges  weakly  to  Y(Sm,k)  where  O'(Sm,i0}jtLi  1S  a  random  walk  on  the  lattice  Gm. 

Proof:  We  have  already  shown  that  T%lk  is  distributed  as  the  n'-m  generation  of  a  branching  process  with 
offspring  generating  function  f(u).  Applying  the  theorem  about  branching  processes  cited  at  the  beginning 
of  this  section  shows  that  15  ~n'~mTn\mik  ~D  15-mlV  also.  For  each  n',  and  m  {27i"jt}?Li  is  a  sequence  of 
independent  random  variables.  Thus  {7t"}£L,  is  also  an  independent  sequence. 

For  any  n',  the  sequence  Yni(S'",  k)  is  a  random  walk  on  G,„.  As  for  {V(Smi*)}jlj,  T(Sj[")  €  G„,,k  = 
1,2,...  with  probability  1.  It  is  straightforward,  although  tedious,  to  show  that  for  sites  a,i,...,rJ,  €  Gm, 
{w  :  y(5"‘,w)  =  X{,i  =  1,. ..,»/)  is  a  continuity  set  for  the  distribution  of  Yt.  The  second  statement  of  the 
proposition  follows. 
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A  consequence  of  this  proposition  is 


Corollary  2.7.  The  sample  paths  o(Yt  ate  uniformly  continuous  in  probability. 

Proof:  The  corollary  follows  from  Skorokhod’s  lemma,  upon  observing  that  for  v  >  v 

p[l|yu  -  mi  <  2y/2-  3-n|n,o  < « < »]  >  <  u  -  t»|n,o  < « < »]  (2.24) 

=  p[t„Xv)  <  u-  v]  (2-25) 

>  1  -  15“n/(u  -  v)  (2.26) 

where  k(v)  is  the  smallest  k  such  that  S„, jt  >  v.  Inequality  (24)  follows  from  the  fundamental  estimate,  (25) 
uses  the  fact  that  {!rn>t}  is  an  i.  i.  d.  sequence  of  random  variables,  and  (26)  is  Markov’s  inequality. 

To  show  that  Y(  is  a  diffusion,  we  must  show  that  it  has  the  strong  Markov  property.  We  consider  those 
stopping  times  of  Y,  that  are  also  stopping  times  of  the  embedded  random  walks  V(S£),  in  an  appropriate 
sense.  Let  C  be  the  class  of  stopping  times  T  of  Y,  such  that 

i.  For  some  n,  Yr  6  Cn  with  probability  1. 

ii.  If  m  >  n  and  T  =  ^*=1  T,,,.*,  then  N„,  is  a  stopping  time  for  the  random  walk  Y(Sm,t)- 
Clearly  S,,^-  €  C  for  any  n  and  k,  so  the  class  C  is  not  vacuous. 


Proposition  2.8.  Let  T  €  C.  Then  Y,  satisfies  the  strong  Markov  property  with  respect  to  T.  That  is  for 
any  bounded,  measurable  functions  f  :  R<J  — ♦  It,  g  :  R*  —  R,  and  S\...sj,  ti...it, 
E[g{Y(siAT),...,Y(sdAT))f(y(1l+T) . V(/*  +  r))] 

(2.27) 

=  £  ff(l'(*,  AT) . Vfo  A T))E); [f(Y(t ,),...,V(4))  ] 


Proof:  By  the  monotone  class  theorem,  it  suffices  to  prove  the  results  for  g(x  i,...,av)  and  f(y\ ,  ...,»/*•) 
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bounded  and  continuous.  For  n  =  1,2,...  let 


yn,t  =  Y,  YiSnfMSn*  <  1  <  ^+l)  (2-28) 

I 

Obviously,  sup,  |y,  -  Y„tt\  <  •  3-n.  Thus,  yn>,  — r  Yt  uniformly  with  probability  1,  as  n  — »  oo.  By  the 

strong  Markov  property  of  random  walks  on  G„  and  the  independence  of  T„t *, 


E[g(Yn(si  AT),...,  Yn(sd  A  T))/(yn(T  +  h), . . . ,  Yn(T  +  t*))] 
=  E  [*(y» (. si  AT),...,  Yn(sd  A  T))Elr  [f(Yn(T  +  <,) . Yn(T  +  <*))]] 


(2.28) 


Since  T  6  Gi  almost  surely, 


E[f(Yn(T+n), . . . ,  y„(T+tt))|y„,r]  =  XX*[/(y„(ti), . . .  ,y„(t*))]  •  p[y„,r  =  *,].  (2.29) 

p 

The  distribution  of  Ynp  is  fixed,  so,  xp  and  P\Ynj  =  xp]  are  fixed  as  n  varies.  Since  /  and  g  are  bounded 
and  continuous, 


<7(y,.(si  A  T), ...,  Yn(sd  A  T))  -  <7(y(s,  A  T), ...,  Y(sd  A  T))  a.s. 

/(y«(T + /,), ....  y„(r + t*))  -  /(y (t + to,  ....  v(r + <*)} 
/(yn(tO,...,y„(to  -*/(y(ti),...,y(<*))  a.s. 


As  ti  — 1  co,  the  dominated  convergence  theorem  gives 


E(g(Y„(s1AT),...,Y„(sdAT))f(Yn(ii  +  T) . Yn(tk+T))} 


-  £fo(y(«,  a  r), ...,  y(sd  a  T))/(y  (tj  +T),...,y(t*  +  t))] 

and 

E\g(Yn(siAT),  ...,Y„(sd  A  T))EYT[f(Yn(U),  ...,!'„(#*))]] 

-  ^[</(y  (si  a  t) . y(s<f  a  t))e$  [/(y  (/, ), ...,  y  (/*))]] 

which  proves  (he  proposition. 


(2.29) 

(2.30) 

(2.31) 


(2.32) 


(2.33) 


We  next  use  this  limited  form  of  the  strong  Markov  property  to  show  that  the  laws  Pr,  x  6  Too  are  uniformly 
weakly  continuous  in  x.  We  do  this  by  studying  escape  times,  the  time  required  for  the  process  to  reach  Gk, 
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starting  from  Gn  for  n  >  k.  Clearly,  escape  times  are  stopping  times  in  the  class  C,  so  the  limited  Markov 
property  in  the  preceding  preposition  applies  to  them.  We  begin  by  proving  a  lemma  for  Markov  chains 

Lemma  2.9.  Let  P  be  an  N  x  N  stochastic  matrix.  Let  Q  be  the  space  of  measures  on  {1, Ar},  with  the 
total  variation  norm.  Then,  the  transformation  p  — ►  pP  is  a  strict  contraction  on  Q  if  and  only  if  no  rows 
of  P  are  mutually  singular. 

Proof:  Let  A  and  p  be  two  distinct  measures  in  Q.  Without,  loss  of  generality,  suppose  that  A j  —  pj  > 
0 ,j  =  ,  Aj  -  pj  <  0 ,j  =  M+  1  ,...,N.  Let  v  =  (*j  ~  Pi)  =  £w+i(Pi  ~  •*;)•  As  p  and  A  are 

distinct,  v  >  0.  Let  votj  —  (A j  —  Pj),j  =  1  ,...,M  and  vfy  =  i'~1(pj  —  A j),j  =  M  +  Then, 


d(XP,  pP)  =  1 £  )PL<  (2-34) 

»  i 

.  M  N 

—  2  E  E(Ai  ~  Pi)pJ,i  ~  E  (P*  “  A*)P*,t-  (2.35) 

»  i=i  *=A/+i 

.  A/  A’ 

^jEE^v-  E  (2-36) 

i  i=l  it=.w+] 

A/  AT 

<^EE  (2-37) 

i=l  *=A/+1  i 

<  i/  max  d(Pj,Pt;)  (2.38) 


Inequality  (2  37)  above  follows  because  f(x,y)  =  |x  -  yj  is  convex  in  both  x  and  y\  the  other  relationships 
are  straightforward.  If  maxj.jt d(Pj,Pk)  <  1,  d(pP,\P)  <  u  =  d(p,\).  If  d{Pi,P3)  =  I  for  some  i  and  j, 
then  take  p  —  6;,  A  =  Sj  to  get  d(p,  A)  =  1  =  d(pP,  A P). 
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Let  x  €  Gj.  Let  Yf  denote  the  continuous  process  constructed  at  the  beginning  of  this  section,  starting 
from  x.  For  A-  <  j,  let  Tk  =  inf :  Yf  €  Gk}.  Then  Tk  is  a  stopping  time  such  that  Y(x,Tk)  £  Gk  almost 
surely.  By  the  strong  Markov  property  of  the  embedded  Markov  chains  Y(Tk,n)  {Y{x,Tk),  k  <  j }  form  a 
discrete  time  Markov  chain. 


With  probability  1,  Y(x,Tk)  will  be  one  of  the  four  corners  of  the  square  Sk  on  level  k  enclosing  x.  Number 
these  four  corners  1,2,3  and  4,  starting  with  1  at  the  Northeast  corner  and  proceeding  clockwise  around 
the  square.  If  Y(x,Tk )  =  i,i  =  1,2, 3, 4,  then  the  transition  probabilities  P[y(x,7i_i)  =  j\Y(x,Tk)  =  t] 
will  depend  on  which  of  the  five  level-A*  squares  of  Sk~ j  x  belongs  to.  Number  the  four  outside  squares 
I,  II,  III,  and  IV,  starting  in  the  Northeast  corner  and  numbering  clockwise  around  the  outer  squares.  Let 
V  denote  the  center  square  of  By  using  the  generating  functions  given  in  display  number,  it  is  not 

hard  to  calculate: 


‘  1  0  0  0  ]  r 1/12  3/4  1/12  1/12] 

p(J)=  3/4  1/12  1/12  1/12  .  (/7)_  0  10  0 

1/2  1/6  1/6  1/6  ’  1/12  3/4  1/12  1/12 

L  3/ 4  1/12  1/12  1/12  J  L  1/6  1/2  1/6  1/6. 

■  1/6  1/6  1/2  1/6  I  r  1/12  1/12  1/12  3/4  1 

!/12  !/12  3/4  1/12  .  p(iv)  _  1/6  1/6  1/6  1/2 

0010’  1/12  1/12  1/12  3/4 

Ll/12  1/12  3/4  1/12 J  0  0  0  1 

■1/2  1/6  1/6  1/6- 
p(\')=  1/6  1/2  1/6  1/6 

1/6  1/6  1/2  1/6 
.1/6  1/6  1/6  1/2. 

Direct  calculation  shows  that 


sup  d(P°(i),  Pa(j))  =  iQ  =  7,  II,  III,  IV  and  sup  d(P^\i),P^(j))  =  I  (2.39) 


Let  BC( r)  denote  the  bounded,  real-valued  continuous  functions  on  T.  For  g  €  BC{T).  x  £T,t  >  0,  let 
P:9(x)  =  E'g(Yt). 


16 


Theorem  2.10.  For  all  t  >  0,  Pt  :  BC(T)  — +  BC(T). 


Proof:  Fix  l  >  0.  For  g  6  BC(r),  we  must  show  that  Ptg(-)  €  BC(T).  To  do  this,  it  will  suffice  to  show 
that  Ptg(x),x  €  GTC  is  uniformly  continuous  and  observe  that  a  uniformly  continuous  function  on  GM  has 
a  unique  uniformly  continuous  extension  to  F. 

Let  x,y  €  Goo,  choose  M  sufficiently  large  that  x,y  €  Gm,  and  suppose  x  and  y  lie  within  the  same  square 
on  level  m  of  Goo,  sav  the  square  Sm.  A  fortiori,  x  and  y  lie  in  the  same  square  Sn  on  level  n  of  Gx.  for 
all  n  <  in. 


Let  e  >  0,  and  let  g  be  a  continuous  function  on  T  with  j|<7||  <  1.  Then 


|£*sQ«)-£W)l 


<  \E*g(Yt)  -  Erg(Yt+T)\  +  |F^0«+r)  -  Eyg(Yt+T) I 
+  \Eyg(Yt)-Eyg(Yt+T))\. 


We  estimate  the  terms  separately. 


Note  that 


EIg(Y(i  +  Tn ))  =  £  Af  W);  Eyg(Y(t  +  Tn ))  =  £  A?  E*g(rn(t)) 


where  Ax  and  Xy  are  the  escape  distributions  of  (/)  on  the  corners  of  S„  starting  from  x  and  y,  respectively. 


Then, 


\ETg{Yn{t  +  n))  -  Ey g{Y„ (t  +  r*))|  <  £  |Af  -  Afli^O'nCO)  <  £  |Af  -  A?| 


Choose  6  >  0  such  that  |u- v|  <  5  implies  |.9(u)-y(r)|  <  j,  and  let  A/,  =  {.<;  :  sup(<u<lI<t+A  IP",, -Vu||  >  6}. 
We  have  shown  that  the  paths  of  Yt  are  uniformly  continuous  in  probability ,  thus,  choose  /i  >  0  so  that 
P[Ah)  <  i-  Then, 


|Fxff(yt)  -  Exg(Yt+T)\  <  £r|<?(V',)  -  y(y+r)| 


<  Ex\g(Yt)  -  s(y,+T)|l([r  >  h)  U  [A%]) 


(2.44) 


+  Ex\g(Yt)  -  g(Yt+T)\l[T  <  t,,Ah]  (2.45) 

<2-|  +  2-|  +  c  (2.46) 

=  3c  (2.47) 

Similarly,  we  can  estimate  | Exg(Y(t))  -  ETg(Y(t  +  Tn))|  <  3c  for  suitable  n. 

To  complete  the  estimate,  choose  n  sufficiently  large  so  that  PX[T„  >  t?]  <  4  Then,  if  we  apply  this  to  the 
preceding  inequalities,  we  get 

\Exg{Yn(t))  -  E»gQrn(t))\  <  3c  +  c  +  3c  =  7c  (2.4S) 

If  we  choose  k  and  m  as  in  the  preceding  paragraphs,  and  let  \x  —  j/|  <  3“m  then  there  must  exist  some 
z  €  Tco  such  that  x  and  z  and  y  and  z  each  lie  within  a  common  square  of  level  m  of  T.  So, 
l*W„(0)  -  &g(Yn(t) |  <  \Exg(Yn(t))  -  Exg(Yn(t)\  +  \E:g(Y„(t))  -  E*g(Yn(t) \ 

(2.49) 

<  14c 

Let  n  — +  oo.  Then  if  |x  —  J/|  <  3_m,  |£x£(y,)  —  Ev g(Y, )|  <  14c. 

This  shows  that  Exg(Yt)  is  a  uniformly  continuous  function  of  x,  for  all  r6rN,  all  /,  and  for  g  an  arbitrary 
uniformly  continuous  function. 

Proposition  2.11.  Yt  is  a  Markov  process. 


Proof:  Fix  t  >  0  and  consider  where  [x]  denotes  the  greatest  integer  less  than  x.  Then,  using  our 


earlier  notation, 


s,  Vf  I15''1!  (,-■  1  '‘f'V  .) 

",15  '  ]5r'  I  []o"(]  P  "Jj 
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As  n  —  oe, 


[15"  <] 

15"  •  - r  V  T„  j  — *  1  a.s.; 

15"t]  f-f  J 

1  J  j=i 


[15"*] 

15n 


—  t 


(2.51) 


The  first  limit  follows  by  applying  the  law  of  large  numbers  to  the  i.  i.  d.  sequence  Tnj-  The  second  is 
simple  analysis.  Thus,  5ni[isnt]  — » i  almost  surely,  asn-*oo. 

If  g(xi,...,Xd)  and  f(yi,—,yk)  are  bounded  continuous  functions  and  sj  <  ...  <  Sd  <  t  <  t\  <  ...  <  tk  then, 
by  Proposition  1.8,  for  any  n, 

E\g(Y(si  A  Sni[i5n<]), . . .  ,y(5d  A  S„t|iB»ti))/(y(fi  +  Sn,[i5»«])> . .  ■ ,  Y(tk  +  •Sn.ps-i)))] 

(2.52) 

=  ^(7(s,  A Sn>[ ,5.|]), •  ■  -,Y(sd  A  Sn,[15Ml))£[/(y(ti),  •  ■  ■  ,y(<t))]|y(5„jis-il)j 
Let  n-*cc,  apply  the  continuity  of  /  and  g,  Theorem  1.10,  and  the  sample  continuity  of  Yt.  We  get 

(2-53) 

This  proves  the  proposition. 


To  show  that  Yt  is  a  strong  Markov  process,  we  will  apply  the  following  theorem: 


Theorem  2.12.  Let  Xt  be  a  process  satisfying 

i.  For  all  f  €  BC(T),  and  allt  >  0,  0  <  /  <  1  implies  0  <  Ef(  X,)  <  1. 
ft.  For  all  s,1  >  0,  Ef( Xt+t)  =  £*(£/(. Xt)). 

Hi.  For  all  f  €  BC(T),  ETf(Xt)  —  /(*)—*  0  uniformly  in  x  as  t  —  0.  Let  7?  =  a  (A',  :  s  <  i)  and 
—  n, > j 7£ .  Then  ifT  is  an  7 . s.  stopping  t/me,  E  rj\7'j'  —  E^  ^  [?/] ,  for  any  measure  p  and  any 

bounded  measurable  y. 

Proof:  This  theorem  is  given  in  Williams  [17],  as  Theorem  III. 15.3. 

To  apply  the  theorem,  it  remains  to  show  that  ETf(Y ()  -  f(x)  —  0  uniformly  in  x  as  t  —  0  for  any  bounded 
continuous  /. 
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Consider  the  operator  Pt  :  BC(T )  — ►  BCIT^)  defined  by  Ptg(x)  =  Ezg(Yt),x  €  Too-  Since  Ezg(Yt )  is 
uniformly  continuous  for  x  £  rTO,  we  can  extend  Ezg(Yt)  to  a  uniformly  continuous  function  defined  for  all 
x  £  Too  and  all  t  >  0.  Let  Ptg  denote  this  extension  of  Ptg.  To  apply  the  theorem  from  Williams  [17]  stated 
earlier,  we  need  only  prove  the  following 

Proposition  2.13.  For  all  f  £  BC(T),  and  all  t  >  0,0  <  f  <  1  implies  0  <  Ptf  <  1,  and  Ptf(x)  —  f(x)  —>■  0 
uniformly  in  x  as  t  —?  0. 

Proof:  Let  /  £  BC(T).  If  0  <  /  <  1,  then  0  <  Ezf(Yt)  <  1  (a.  s.),  since  conditional  expectation  is  a 
positive  operator.  Since  I\  BC(r),  and  [0,  oo)  are  all  separable,  we  can  modify  P  on  a  single  null  set  N  to 
get  0  <  Ezf(Yt)  <  1. 

We  apply  the  following  lemma. 

Lemma  2.14.  Let  Pt  :  BC(T)  —  BC(F)  be  a  substocbastic  Markov  semigroup.  If  Ptf(x)  —  f(x)  as  t  J  0 
for  all  x  £  F  and  f  €  BC(T),  then  Ptf(x)  -  f(x)  — *  0  unifornily  in  x  as  t  — ►  0,  for  all  f  £  BC(T). 

Proof:  This  paraphrases  formula  III.S.2.iv.  in  Williams  [17]. 

Let  e  >  0.  Let  x  £  rTO,  and  let  g  £  BC(r).  Given  e  >  0,  choose  5  such  that  |y  -  x\  <  6  implies  that 
|<7(r)  -  ^(y)|  <  4-  Fr[|y<  -  x|  >  5]  j  0  as  t  J  0.  Choose  t  sufficiently  small  that  Pz[\Yt  -  r|  >  6]  <  Then 

l-P:S(^)-5(x)l 

<  1^(50:)  -  9(x))W  ~  x\  <  6) |  +  |£*(p(yf)  -  5(x))l[|y<  -  x\  >  6)\ 

(2.54) 

<  iPT[\Yt  -  x\  <  f\  +  2  •  Fr[|y«  -  x\  >  6) 

<f+f=£ 

This  completes  the  proof  of  the  proposition. 
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In  terms  of  Williams  [17],  this  makes  Pt  an  FD  semigroup  and  Yt  an  FD  process.  Applying  the  theorem  that 
we  stated  earlier  in  this  section  to  Yt,  we  conclude 

Corollary  2.15.  Yt  has  the  strong  Markov  property. 

Henceforth,  we  will  refer  to  Yt  as  the  snowflake  diffusion  . 

In  fact  we  can  show  that  the  snowflake  diffusion  is  the  unique  diffusion  limit  for  any  of  the  random  walks 
Xp,  within  a  constant  change  of  time  scale 

By  Proposition  2.6,  random  walks  on  all  the  lattices  Gn  are  embedded  in  the  snowflake  diffusion.  We  also 
note  in  passing  that  for  Vo  =  0,  3 Yt  —n  Fist;  we  will  return  to  this  point  in  the  next  section  at  greater 
length.  These  two  facts  form  the  basis  for  the  proof  of  the  following  proposition. 

Proposition  2.16.  Within  a  non-random  change  in  time  scale,  Yt  is  the  unique  limit  in  distribution  of  the 
processes  Yn(t). 

Proof:  For  each  n,  define  the  sequence  T” ,  7? , . . .  as  in  proposition  2.6.  for  each  k,  Tf  is  the  sum  of  k 
independent  random  variables  equal  to  15~nHr  in  distribution.  If  Var(W)  =  o'2,  then  ETf  =  k  •  15-'1  and 
Var  {If)  =  ko2. 

For  arbitrary  n  and  t, 

Yn{t)  =  3-"A'([l5"f])  =d  3-y(r  ([15°0))  =  Y  (lo-nr’([15"t])) .  (2.55) 

Note  that  t  —  15-n  <  15-n  -  [I5"f]  <  t  +  15“".  Now,  T£  is  the  sum  of  k  independent  random  variables,  with 
distribution  15-m  ■  W.  Thus, 

E(lo~'n  ■  T([15m/]))  =  15-m  ■  £(r([l5m<]))  =  15-m  ■  [15‘m  (2.56) 
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<r(15-m  •  !T([15mt]))  =  lo~~m  ■  <r(T([15mi]))  =  lo~2m  ■  15"m  •  <r2  =  15“m  • 


0  (2.57) 


as  n  — ►  oo,  where  <r2  =  <r~(W).  Thus,  as  n  — *■  oo,  Y„(i)  ~~p  Yt. 

By  the  same  argument,  for  any  fi  <  ...  <  Ik  and  any  convergent  subsequence  Y„»{t)  converges  in  distribution. 
Then  by  the  same  analysis  as  in  the  preceding  paragraph, 

-+D  (Y(U),...,Y(tm))  (2.58) 


Thus  all  limits  of  {yn(t)}^j  have  the  same  finite  dimensional  distributions.  Since  the  finite  dimensional 
distributions  of  a  stochastic  process  determine  the  process,  it  follows  that  Yt  is  the  unique  limit  of  the 
processes  {y„(t)}~=1. 

Now  we  consider  again  the  general  model  for  random  walk  proposed  at  the  beginning  of  this  section.  Let 
Xp  be  random  walk  on  the  snowflake  lattice  with  parameter  p.  Let.  po  =  p  and  for  n  =  1,2, . . .,  let  p„  be  the 

probability  that  Xp  starts  from  0  and  reaches  (3”,3n)  before  either  (3",0)  or  (0,3").  It  is  not  hard  to  see 

that  in  general  pu  pn-\-  We  can  find  a  recurrence  relationship  by  calculating  computing  p\  as  a  function 
of  po  Observe  that  we  can  identify  vertices  of  U\  that  are  reflections  of  one  another  across  diagonals  of  the 
square.  Thus,  we  can  identify  the  following  sets  of  vertices 

a  =  {(0,0)}  6  =  {(0,1), (1,0)}  c=  {(1,1)} 

d=  {(1,2),  (2,1)}  e  =  {(0,2),  (1,3),  (2,0),  (3, 1)}  /  =  {(2,2)}  (2.59) 

9  =  {(2,3),  (3,2)} 

For  each  i  €  {a,  ■ . .  ,<?},  let  g,  be  the  probability  that  X„  starts  from  \ertex  i  and  reaches  (3, 3)  before  either 
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(0,3)  or  (3,0).  Using  the  Markov  property  of  Xr,  we  get  the  following  system  of  equations  for  q, 
qa  =  (l  -  p)?t  +  pqc  qt>  =  5(1  -  p)qa  +  pqb  +  k  (1  -  p)qc 

qe  =  \pqa  +  5(1  -  p)qb  +  |(i  -  p)qd  Jr  ^pqj  q<i  =  K1  -  P )?« +  hw*  +  K1  ~  +  *  (*  -  p)<u 


<?e  =  |(1  -  P)qd  +  jWe 


9/  =  \pqc  +  5(1  -  p)9<j  +  -K1  -  p)?s  +  \p 


qg  =  K1  -  P)«/  +  P?s  +  K1  -  p) 


Applying  Cramer’s  rule,  we  get 


(4)5-(l-p)3-(l  +  p)3  1 

(|)5  ■  (1  -  P)3  •  (1  +  P)3  •  (4  -  3p)  4  -  3p 


Formula  (2.61)  gives  us  a  recurrrence  for  pi  in  terms  of  po-  Furthermore,  we  see  that  p  =  5  is  a  fixed  point 
for  this  recurrence  formula,  and  that  for  any  po  <  1,  pn  — ►  5  as  n  — i ►  00. 

For  73  =  1, 2, . . .,  let  Tn,p  =  T(n,p)  be  the  number  of  steps  Xp  requires  to  pass  from  0  to  some  other  vertex  in 
3"  •  G.  In  analogy  with  the  branching  property  of  A'1/3,  it  is  not  hard  to  see  that  Tn,r  =  T£~l'p,  where 

T£~1,pk  =  1,2,...  is  a  sequence  of  random  variables  distributed  as  rR-1,f>,  B  is  distributed  as  r(l,p„_i) 
and  B  and  T£~l’pk  =  1,2, . .  .  are  mutually  independent. 

Let.  m„  =  £T(l,p„),  and  let  M„  =  f]"=1  m(-  Then,  ETn,p  =  Mn,  and  an  argument  analogous  to  that  for 
branching  processes  shows  that  Tn,p  ■  M~l  is  a  positive  martingale  converging  a.  s.  and  in  distribution  to 
some  random  variable  W. 


Lemma  2.17.  Within  a  constant  scaling  factor,  IV  =#  IV. 


Proof:  Let 


/»(«)  =  EuT^'p-\  gn(\)  =  Eex p  (~XT(n,p)  ■  M~') 


Then  from  the  branching  property  of  Tn,p  outlined  in  the  preceding  two  paragraphs 


gn(X)  =  fn  (ffn-i(X-r n,T 1 )) 


(2.63) 
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Since  Tn'p  ■  Mn  1  —*d  W,  gn( A)  — -  #(A)  =  jE'exp(-Aiy)  as  n  — ►  oo.  pn  —■ ►  1/3,  so  fn( u)  — *■  /(u),  and 
i/i„  —  15.  Substituting  these  limits  in  (2.63)  gives 

0(A)  =  /(0(A  •  15-1)),  (2.64) 

the  same  functional  equation  satisfied  by  W .  It  is  well  known  that  the  solution  to  Abel’s  functional  equation 
is  unique,  except  for  a  constant  multiplier  of  the  argument.  (  See,  for  example,  Seneta  [16],  Theorem  3.1.) 
This  implies  that  W  =d  W  within  a  constant  change  of  scale. 

Let  yP(<)  =  3_nA'p([A/nt]).  Arguments  analogous  to  those  in  Theorem  2.2  ,  Lemma  2.3,  Propositions  2.6, 
2.8,  2.11,  and  2.13  show  that  we  can  find  a  diffusion  Y*  on  T  and  a  subsequence  n'  such  that  Yp  —  Y  weakly. 
Let  { T„ }  be  the  sequence  of  times  when  Yt  visits  distinct  vertices  of  C\.  Note  that  T  =d  H'.  As  pn  — »  1/3, 
{Y(Tn)}  =D  {A'1/3} 

Theorem  2.18.  Yt  =d  Yt,  within  a  constant  change  of  time  scale. 

Proof:  As  a  consequence  of  Lemma  2.17,  we  can  choose  a  time  scale  for  VJ  so  that  3Y«  =d  Yist-  The  proof 
now  proceeds  as  in  Proposition  2.16.  Let 

y*(0  =  3-my(T([15m/]))  =  y(15-m  ■  f([15m/]).  (2.65) 

As  Yt  has  continuous  sample  paths,  1~(<)  — r  Yt  almost  surely.  On  the  other  hand,  }'„’(/)  Ym(i),  and 
y,j)(  1)  —d  Yt ■  Thus,  Yt  =D  Yt.  This  proves  the  theorem. 
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3.  Scaling  Properties: 


In  this  chapter,  .ve  will  study  some  of  the  detailed  sample  path  properties  of  the  snowflake  diffusion.  Let 
n  >  0,  let  x  and  y  be  points  on  Gx  D  [0,3]  such  that  x,y  £  S  for  some  square  S  on  level  n.  Let 
Ty  =  inf{t  :  Yt  =  j/}  be  the  first  visit  of  the  diffusion  on  Fj  to  y\  set  Ty  =  oo  if  this  set  is  empty. 

We  begin  by  proving  two  technical  lemmas. 

Lemma  3.1.  ExTy  <  369  ■  3-”. 

Proof:  Suppose  x,y  £  S.  There  exists  m  >  n  such  that  x,y  £  Gm.  We  will  establish  our  estimate  by 
induction  on  m. 

Suppose  that  m  =  n.  Then  x  and  y  are  necessarily  corners  of  5.  Let  T\,  Tj, ...  be  the  times  between  successive 
visits  to  G„  Then,  Ty  =  Tj  where 

:( 

M  is  the  number  of  corners  of  S  that  Yt  visits  before  Ty,  N,  is  the  number  of  excursions  to  Gn  \  S  between 
the  time  Yt  hits  the  i  -  1“  and  itb  distinct  corners  of  S,  and  the  number  of  points  in  G„  \  S  that  Yt 
visits  on  the  j1h  such  excursion. 

The  strong  Markov  property  shows  that  M  has  a  geometric  distribution  with  parameter  For  i  =  1,2,... 
Ni  either  is  identically  0  or  else  has  a  geometric  distribution  with  parameter  To  estimate  ER,t},  note 
that  each  excursion  outside  of  S  is  a  random  walk  on  a  finite  graph,  and  JR,  }  is  the  number  of  steps  the 
walk  takes  to  return  to  its  starting  point.  It  is  well-known  that  this  expected  time  is  equal  to  twice  the  total 
number  of  edges  in  the  graph  divided  by  the  degree  of  the  starting  vertex;  see  Gobel  and  Jagers  [9]  for  this 
and  other  general  results  about  random  walks  on  finite  graphs.  Thus,  ER,j  is  proportional  to  the  number 
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of  (he  edges  in  the  graph  cut  out  of  G„  by  S.  This  number,  in  turn,  is  less  than  the  number  of  edges  in  Gm. 
It  is  not.  hard  to  compute  that  Gn  has  6  •  5n+1  edges,  so,  ERij  <  20  ■  5n 


Recall  that  ETi  =  15-n,  and  apply  Wald’s  identity. 

ExTy  =  EM  •  (. EN(i )  •  ERij  +  1)  •  ET 


<  3  ■  (2  •  20  •  5"  +  1)  •  15-n  (2.2) 

<41 ,3-(n-i) 

V  now  proceed  by  induction.  Suppose  that  if  x,y  6  Gm  with  |x  -  y|  <  3-n,  then  ExTy  <  82  •  ^2™  3”,+1. 
Let  x,y  €  Gm+ j,|z  —  y|  <  3"”.  There  exist  points  v  and  w  in  S  n  such  that  x  and  v  and  y  and  w 
are  adjoining  corners  of  squares  on  Gm+\.  (  x  and  v  may  be  identical,  and  so  forth.  )  The  strong  Markov 


property  gives 


EzTy  <  ETTV  +  EVTW  +  EwTy 


m-1 

<  41  •  3_m  +  82  •  3~*  +  41  •  3-m 

n-l 

m+1 

=  82  -J2  3-'+1 

n 

completing  the  inductive  step. 


(2.3) 


Finally,  if  x,y  €  <S  n-oo,  then 

CO 

ErTy  <82-^3“f+1  =  369-3“n  (2.4) 

n 

This  completes  the  proof  of  the  lemma. 


Corollary  2.2.  If  a-,y  e  n  [0,3],  such  that  |x  -  y|  <  3-'1,  then  ET7'y  <  738  •  3~n. 

Proof:  If  |x  -  y|  <  3~",  then  either  they  lie  in  a  single  square  S  or  level  n,  or  else  the\  lie  in  two  such 
adjoining  squares.  In  either  case,  the  preceding  lemma  gives  the  desired  bound. 

Recall  that  3  G„  =  G„_i  This  fact  and  the  uniqueness  of  the  snowflake  diffusion  limit  imply,  the  following 
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Lemma  3.3.  Let  x  €  T.  Then  3 Yf  =d  Y{ 

Proof:  Let  x  6  Using  the  notation  of  section  2.  consider  the  sequence  of  processes  3  -  Yn(x.i).v  =  1K. 
We  have  shown  that  Y„(x,t)  —  Yf  as  n  —  oo,  so,  trivially.  3  -  Y„(x,l)  —  3  •  V* 

On  the  other  hand, 

3  •  Yn(x:i)  =  3  •  3-nA'(3nx,  [15”t])  =  3-n+1A'(3n-13x5  [15"_1(15/)])  =  V;_i(3x,  15t)  (3.4) 

Thus,  3  -  Yn(x,<)  —  Y^j .  Since  both  of  these  limits  are  unique,  it  follows  that  3}'*  =d  Y^,  for  x  6  Tco- 

For  general  x  in  T  choose  a  sequence  x„  in  Too  converging  to  x.  Then,  for  each  n,  3Y,X"  =&  Y^,"-  As  Y  has 
the  Feller  property. 

3YX"  —d  3Ytx  and  Yt3r"  —D  Y?T  n  —  oo  (3.5) 

Therefore.  3Ytx  —D  Vis*. 

Corollary  3.4.  3Y°  =D  Y°5r 

For  starting  points  in  other  than  0,  a  result  similar  to  the  preceding  corollary  also  holds.  Let  x  6  Tx. 
and  define  M,  :  y  —  3_1  ■  (y  -  x)  +  x.  Then  Mxx  =  x. 

Corollary  3.5.  Let  x  €  Gx.  Then  there  exists  a  stopping  time  Tt  >  0  such  that  MrYtT  —d  Yfbt ,  0  <  i  <Tt 

Proof:  Suppose  that  x  €  G„.  Then  there  exists  either  one  or  two  squares  S  on  level  n.  such  that  x  is  a 
corner  of  S  However,  if  x  is  a  corner  of  5,  then  there  is  some  square  S'  on  level  n  +  1  contained  in  S  such 
that  x  is  also  a  corner  of  S'.  Clearly,  M.  :  S  —  Si.  Let  V  be  the  union  of  any  squares  S'  on  level  n  -f  1 
adjoining  x  Then,  the  proof  of  the  fractal  scaling  law  will  continue  to  work  if  we  set  Tx  to  be  the  first  time 
the  process  started  from  x  exits  L"  and  Tx  satisfies  the  requirements  given  in  the  statement  of  the  corollarv. 
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We  shall  call  the  scaling  property  stated  in  either  corollary  the  unbounded  /racial  scaling  law. 


If  we  restrict  our  attention  to  Tj.  then  our  scaling  properties  are  not  quite  so  tidy.  However,  a  closely  related 
property  does  apply,  which  arises  from  the  invariance  of  T*  under  a  family  of  contractions.  Define 
Ari  :  x  — r  3  •  x 


jVo  :  x  — —  3  •  x  +  (6, 6) 
Ar3  :  i  —  3  •  x  —  (6, 6) 

N4  :x— 3-*-(U,4§) 


Ar5  :x-‘3-*-(4§,  if) 


and  Ar(x)  =  < 


r-viM. x  €  [o,  i]2 

N2(x),  x  G  [1 . 2]2 
N3(x),  x  €  [2, 3]2 
A%(*),  x  6  [0, 1]  x  [2,3] 
l  Ns(x),  x  G  [2. 3]  x  [0, 1] 


(3.6) 


Note  that  N  is  a  continuous  mapping  of  Tj  onto  r& 


Lemma.  Let  igT.  Then,  N(Yf)  =D  Y$x\ 

Proof:  Inspection  will  verify  that  if  A'n  is  a  random  walk  on  starting  from  x,  then  AT(A*)  is  a  random 
walk  on  G?*_i  starting  from  Ar(x),  for  any  kl.  We  have  shown  in  Section  2  that  Yt  weakly.  Then 

rnx{u.t])  =  A'(A-,*S.., -  y,5,  (3.T) 

and  since  N  is  continuous,  Ar(A'^sV<j)  —*  N(Yt). 

We  shall  call  this  scaling  property  the  bounded  /racial  scaling  law.  As  an  important  consequence,  note  that 
for  any  x  g  Fj,  and  any  measurable  set  A,  PJV(r)[y15,  g  A]  =  PT[Yt  €  A'-1  (.4)). 

Our  next  step  is  to  show  that  two  independent  copies  of  the  snowflake  diffusion  restricted  to  Tj  meet  in  finite 
time  with  probability  1.  First  we  prove  two  technical  lemmas.  For  x.y  g  T.  let  d(x.y)  denote  the  Euclidean 
distance  between  x  and  y.  Let  D[ 0,  oc]  denote  the  space  of  functions  u>  .  R+  —  T  which  are  right  continuous 
and  have  left  limits  for  all  /  >  0. 
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Lemma  3.7.  For  any  t ,  c  >  0  and  any  compact  set  K.  let 


H(i,c )  =  {«  :  inf  d(ui{v),K)  >  c,0  <  v  <  t  ~  <}. 

Then  H(i,  c)  is  an  open  subset  in  the  topology  T  on  D\ 0.  oc)  defined  by  convergence  in  the  Skorokhod  metric 
on  compact  internals  [O.p]. 

Proof:  Let  u.v  G  L>[0,co],  For  any  p  >  0,  let  pP(p.  v)  equal  the  infimum  of  those  t  >  0  for  which  there 
exists  a  continuous,  increasing  function  A  :  [0,p]  —  [O.p]  such  that 
i.  sup{0  <i<p:  |A(t)  —  t  |}  <  c 
ii.  sup{0  <i  <p:  d(p(A(t)),J/(0)}  <  c 

It  is  easy  to  see  that  pP  is  a  pseudo-metric  and  that  the  topology  T  is  indue.  b\  the  famih  of  pseudo-metrics 

Suppose  u  G  H(t,c).  Let  a  =  inf{d(«(u),  K ),  0  <  u  <  i  -f  e}  By  hypothesis,  a  —  c  —  6  >  0.  Choose 
v  €  jD[0,oo]  with  pr(w,r/)  <  There  exists  a  strictly  increasing  continuous  function  A  ;  [0,<]  —  [0,/] 
satistfying  i.  and  ii.  By  the  triangle  inequality,  for  any  v  G  [0,<], 

d(v(u),  K)  >  d(u{ A(«)),  A')  -  </(*( A(u)),  */(«))  >  e  +  6  -  |  >  e:  (3.S) 

since  A  :  [0,t]  — » [0,/],  d(«(A(u)),  K)  >  e  +  6,  for  0  <  u  <  t.  Thus,  v  G  H(i.  ().  so  Ht  is  open. 

Let  T  =  inf{t  :  u(t)  G  A'}.  For  every  e  >  0,  H(l,c)  C  {*>  :  T(u)  >  /}.  Let  C[0. oc]  denote  the  continuous 
functions  from  R+  to  Ti. 

Lemma  3.8.  {w  G  C[0,oo] :  T(u)  >  f}  =  Uf=lH{i,k~l) 

Proof:  Suppose  us  is  continuous,  and  that  for  every  k  >  0,  there  exists  s*  G  [0,  /]  such  that  d(  Us 

As  [0,  <]  is  compact,  there  exists  a  convergent  subsequence  $**  with  s*>  —  sG  [0,#].  As  u>  is  continuous. 
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~{xc>)  —  ^-(x).  However,  d(a;(xi/).A')  —  0,  so  u(s)  €  K,  and  T(n)  <  s  <  i.  Thus,  if  >  t  then 
inf {d(^(s).  A'),  0  <  s  <  t}  >  it-1,  for  some  k. 

Using  the  continuity  of  w  again,  there  must  be  some  m  such  that  d(w(s).  A')  >  <  s  <  /  -j-  m-1.  If 

not.  there  would  exist  a  sequence  x„  —  t  such  that  dfe(xv).K)  <  k~l.  forcing  d{ps{i),K)  <  krl.  Choose 
N  =  min{Jt,  m}  to  give  w  6  //.y-j  - 

To  restate  the  previous  result 

Corollary  3.9.  T :  C{9,oo)  —  [0,oc)  is  a  lower  semicontinuous  function. 

Proof:  Since  Ufiji/i-i  is  an  open  set,  this  follows  by  definition. 

We  now  proceed  to  prove  the  main  result  of  this  section. 

Theorem  3.X0.  Let  Yt  and  Y{  be  two  independent  copies  of  the  diffusion  on  IV  Let  Tm  =  inf{u  :  Yv  =  VJ}. 
Then  Tm  <  oo  a.  s. 

Proof:  In  Aldous  [l]  it  is  shown  that  if  A",  and  A','  are  two  independent  copies  of  a  continuous  time  random 
walk  on  a  finite  graph  H,  then  there  exists  some  constant  D  such  that  ETu  <  Z?max.j  E'Tj.  where 
the  maximum  is  taken  over  all  pairs  of  states  i,j.  By  Markov’s  inequality,  it  follows  that  P[T\{  >  /.]  < 
i~l  ■  Dmax;jE'Tj. 

Let  VB(/)  be  a  random  walk  in  continuous  time  on  Cn-  We  have  shown  in  the  preceding  section  that  for  any 
pair  of  vertices  x  and  y,  such  that  j|x  —  y||  <  3_I  then  E'T}  <  738  ■  3”’.  If  we  regard  }'„(/)  as  a  sequence  of 
processes  on  Tj,  then  P[T^  >  t]  <  t~l  ■  13SD  for  all  n. 

Consider  the  set  A  =  {(x,x)  :  x  €  I^}.  Clearly,  Tu  is  the  first  hitting  time  on  A  for  the  process  ()'»,}’/). 
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As  before,  let  Ht  =  (w  :  inf  d{y:[v),  A)  >  c.O  <u  <i -i- c.  Since  Ht  C  {X\/  <  #}.  it  follows  that 

Pp'n(0  6  H<)  <  P[Ty,  >i)<  r 1  -  73SD  (3.9) 

for  all  n.  As  //£  is  open,  Prokhorov’s  theorem  and  the  Lemma  3.1  show  that 

P[Y,  €  Ht]  <  liminfP[rn(/)  €  H<]  <  t~'  -  73SD  (3.9) 

n — cc 


As  Yt  has  a.  s.  continuous  sample  paths,  the  preceding  equation  and  Lemma  3.S  show  that 

P[Tm  >  tj  =  lim  inf  P[y  6  //„-.}  <  t~l  -  738 D  (3.10) 

n— -co 

Letting  i  —  oc  gives  P[Tu  =  oo]  =  0,  which  is  what  we  proposed  to  prove. 

We  next  consider  invariant  measures  for  Yt.  For  the  moment,  restrict  the  random  walks  Y„(i)  to  IV  (Or, 
equivalently,  restrict  them  to  [0,3]2).  Since  each  Y„  is  now  a  random  walk  on  a  finite  graph,  it  has  a  unique 
reversible  stationary  measure,  which  we  shall  call  //„.  Regard  yn  as  a  revalued  stochastic  process,  and  //„ 
as  a  measure  on  IV  Since  Tj,  is  compact,  {//„}  has  a  weakly  convergent  subsequence 

Proposition.  Ys  has  a  stationary  distribution  p,  and  p„>  — p  P- 

Proof:  To  show  the  proposition,  we  apply  weak  convergence.  We  have  previously  shown  that 

£[r»(«)  -  y»(s)  ip  r„(o  -  v;(u)ip)  <  d(i  -  s r  (3.m 

for  any  p, -•  >0  and  any  s  <  u  <1.  Since  Tj  is  compact,  the  sequence  //„  is  weakly  precompact.  Let  n'  be 
a  sequence  of  integers  such  that  //„<  converges,  say  to  //.  If  we  let  Y£  be  the  stationary  version  of  Y„>  then 
standard  results  on  weak  convergence  (See,  again,  Theorem  15.6  in  Billingsley  )  show  that  converges 
weakly  to  a  process  Y^,  where  Yl‘  is  a  version  of  Y  with  stationary  distribution  p. 
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Theorem  3.11.  The  distribution  ofY-  converges  p  in  total  variation  norm 


Proof:  Let  p:(y..4)  =  Ps[ Yt  €  A).  Let  Y(  be  an  independent  stationary  diffusion  on  T&  Couple  Yt  to  YJ 
by  letting  y  and  Y'  move  independently  prior  to  Tv  hut  specifying  that  they  move  identically  afterwards. 
Then  for  any  measurable  .4, 

M-4)  -  p(.4)|  <  P[Yt  t  Y'\  <  r1  -  738 D  (3.12) 

and  the  inequality  is  uniform  over  all  measurable  .4.  The  theorem  follows. 

Corollary  3.12.  fi  is  the  unique  stationary  distribution  for  Yt. 

Theorem  3.13.  /i  is  normalized  Hausdorff  log$  3-dimensional  measure,  restricted  to  Ij. 

Proof:  Since  p  is  a  stationary  measure,  we  apply  the  preceding  theorem  and  the  bounded  fractal  scaling 
law  to  get 

f%:  G  Ar-1(4)]  -  pCjV-1^))  (3.13) 

PQ[V,  G  A'"1  (4)]  =  J%15,  g.4]_  p(4)  (3.14) 

as  t  —  oo.  So,  p  satisfies  the  equation  p(4)  =  p(Ar-1(.4)).  Theorem  4.4.1  in  Hutchinson  [12], shows  that 
there  is  a  unique  measure  on  ft  that  satisfies  this  equation.  Since  p-dimensional  HausdorfT  measure  restricted 
to  Tj  also  satisfies  this  equation,  it  follows  that  p  is  p-dimensional  Hausdorff  measure  on  Tj. 

Now  consider  Yt  on  F.  If  we  let  p  be  Hausdorff  log3  5-dimensional  measure  restricted  to  T,  we  have  the 
following 

Corollary  3.14.  p  is  an  invariant  measure  for  Yt . 

Proof:  Normalize  p  so  that  =  1,  and  let  /  be  a  continuous  function  on  T  with  compact  support- 
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Choose  N  sufficiently  large  so  that  T  n  [0, 3A  ]  contains  the  support  of  /.  Then  for  n  >  N,  //  is  an  invariant 
measure  for  y/n)  .  the  diffusion  process  restricted  to  m[0,3*v].  Thus,  as  n  —  co, 

Jfdn  =  J  EW<n)Mdx)  -  J  E^UO'tMdz)  (3.15) 

As  the  continuous  functions  with  compact  support  are  dense  in  L1,  the  result  follows. 

We  have  shown  in  Theorem  3.11  that  converges  io  Hausdorff  log3  5-dimensional  measure  in  total  variation 
norm  for  restricted  to  Ts-  This,  together  with  the  scaling  laws,  implies  that  Y:  has  a  transition  density 
with  respect  to  Hausdorff  measure. 

Recall  the  transformations  A71; ...,  Ar5  and  N  defined  in  (3.6),  Let  /it(  )  =  P-\Y<  €  -}  and  let  //(  )  denote 
Hausdorff  log3  5-dimensional  measure,  restricted  to  T. 

Theorem  3.15.  //,  -C  ft  for  all  i  >  0. 

Proof:  We  begin  by  proving  the  theorem  for  Yt  restricted  to  IV  Let  B  C  Tj  be  a  set  with  p(B)  =  0. 
and  suppose  that  fit(B)  =  q  >  0.  We  observe  in  passing  that  fi(B)  =  0  iff  p(AT(B))  =  0.  and  that 
Ar-1(.V(S))  D  B.  Let  =  |Jf  A7*(P).  Since  p(B)  =  0,  /i(Nk(B))  =  0  for  all  k,  so  p{Bx)  =  0.  On  the 
other  hand,  B  C  A r~k(Bos)  for  all  k.  By  the  bounded  fractal  scaling  law, 

=  /tj(A  (Poo))  ^  /h(P)  =  9  (3.16) 

Thus.  liminfr_c^//j5i,(Pco)  >  q.  But  this  contradicts  the  fact  that  //,  —  /;  in  total  \ariation  norm  as 
1  —  co.  By  contradiction,  fi:(B)  =  0.  Thus,  <  //. 

To  sliow  absolute  continuity  for  the  unbounded  process,  let  Yt<n'  denote  the  diffusion  process  restricted  to 
Fu  [0,3"].  As  n  —  co,  V-/"1  —  in  total  variation  norm.  The  proof  for  V«  restricted  to  Ft  shows  that 


€  ]  is  absolutely  continuous  with  respect  to  HausdorfT  measure  for  all  n.  Therefore,  it  follows  that 


* 


P[Yt(n) 

P[Yt  6  ]  is  also  absolutely  continuous  with  respect  to  HausdorfT  log5  3-dimensional  measure  on  T. 


3-1 


4.  Computing  Generating  Functions 


We  can  determine  the  generating  function  of  the  distribution  of  T 1  by  elementary  calculations.  Let  k(u)  = 
E°vr  where  r  is  the  time  required  to  cross  from  Oto  (1,1).  r  has  a  geometric  distribution  with  p  =  1/3.  so 
k(u)  =  u/(3  —  2u),  trivially. 

Suppose  A'o  =  0,  and  let  a  denote  the  first  hitting  time  on  an  outer  corner,  other  than  0.  Let  /(u)  =  E°ua . 
We  calculate  the  generating  functions  of  some  hitting  times  for  corner  1,  as  a  preliminary  to  calculating  the 
generating  function  of  a. 

Let  p  denote  the  first  hitting  time  on  1,  on  the  on  the  set  where  lis  the  first  outer  corner  of  l'o  that  the 
walk  visits.  As  the  distribution  of  the  random  walk  is  not  affected  by  symmetries  of  l'o,  the  g  depends  only 
on  the  graph  distance  between  the  random  walk’s  starting  point  and  1.  Thus  we  can  identify  vertices  of  Tj, 


a  =  {(0,1), (1,0), (0.2), (1.3), (2,0). (3,1)} 
b=  {(1,1),  (1,2),  (2.1)} 

c={(2,2)} 


d={(2,3),(3,2)} 

For  i  a  vertex  in  U0.  let  <?,-(«)  =  E'v1'.  g  then  satisfies  the  following  system  of  equations 


9a  =  %V9o  +  im 


9b  =  lu9a  +  +  %vgc 


9c  =  \u9b  +  \ugd  +  ju  gd  =  jugc  +  \ugd  +  £u 
Direct  calculation  shows  that  this  system  of  equations  has  the  solutions: 


3(2  —  u)(18  —  iou  +  u2) 
2u(3  —  2u) 

(2  —  u)(18  -  15u  4  u2) 


tr( 3  -  u) 

3(2  -  u)(18  -  15u  4  u2) 
t/(12  -  10»  4  u2) 

(3  —  u)(18  —  low  4  u2) 


Now,  let  A  be  the  first  hitting  time  for  anj  comer  other  than  0.  Again  we  can  use  symmetry  to  identify 
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vertices,  in  the  groups 


^  =  {(0, 1),  (1,0)> 

*  =  {(1,1), (1,2), (2,1)} 
y=  {(1.2),  (2,1),  (2, 2)} 

*  =  {(0, 2),  (1,3),  (2, 0),  (3, 1),  (2, 3),  (3, 2)} 


(4.4) 


Then  /t,(u)  =  E'vx  satisfies 


with  solutions 


kw  =  5 vhw  +  ^uhx  hx  -  \uhx  +  5 uhy 

hy  =  \uhw  +  5 vhy  +  %uh.  +  h2  =  ^uhy  +  | uh;  +  5U 

u2(3  —  u ) 


(4.5) 


(2  -  u)(18  -  15u  +  u2) 

L  _  2u(9  —  3 u  —  tr) 

— 


hx  = 
h2  = 


(2  -  u)(18  —  15u  +  u2) 
u(36  —  30u  +  5u2) 


(4.6) 


v  3(2  —  u)(18  -  15u  +  u2)  "z  3(2-u)(18-  lou-fu2) 

Using  the  generating  functions  g,  h,  and  k ,  we  can  now  calculate  /(«).  Start  at  Oand  condition  on  N,  the 
number  of  returns  to  Obefore  hitting  a  different  outer  corner.  We  use  the  symmetry  of  Vo  and  the  generating 
functions  we  have  just  calculated,  to  get 

M») 


=  *(«)£Mu)MtO*(u)]"  = 


(4.7) 


If  we  substitute  the  generating  functions  we  have  already  calculated,  we  get 

|  _ T3 

VI  /TV  W-  TV 

u 


Euua  = 


3  —  2u 


1  - 


(2  -  u)(18  —  15u  +  u2) 
2«(3  -  2u) 


i/ 


(2-u)(18-15u  +  u2)  (3-2«)J 


(3  -  2t/)(12  -  12u  +  u-) 

=  /(«) 


(4.8) 

(4.6) 

(4.9) 


The  symmetry  of  Uq  shows  that  f(u)  is  the  generating  function  for  the  distribution  of  the  time  to  cross 


between  any  two  distinct  corners  of  G. 
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